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Abstract. Starting from the semi-classical spectrum of Schrodinger 
operators — h 2 A + V (on R n or on a Riemannian manifold) it is possible 
to detect critical levels of the potential V. Via micro-local methods one 
can express spectral statistics in terms of different invariants: 

• Geometry of energy surfaces (heat invariant like). 

• Classical orbits (wave invariants). 

• But also classical equilibria (new wave invariants). 

Any critical point of V with zero momentum is an equilibrium of the 
flow and generates many singularities in the semi-classical distribution 
of eigenvalues. Via sharp spectral estimates, this phenomena indicates 
the presence of a critical energy level and the information contained 
in this singularity allows to reconstruct partially the local shape of V. 
Several generalizations of this approach are also proposed. 

Keywords : Spectral analysis, P.D.E., Micro-local analysis; Schrodinger 
operators; Inverse spectral problems. 



1. Introduction. 

1.1. Background and basic definitions. In this article we are here in- 
terested in the inverse spectral problem for partial differential operators and 
pseudo-differential operators in the semi-classical or high-energy regime. A 
natural question is to try to understand how the semi-classical spectrum 
of such an operator can describe the shape of the graph of the (principal) 
symbol: critical points, extrema and associated local Taylor expansions. 
Because the spectrum is invariant under translation of the symbol it is in 
general not possible to obtain more than a qualitative answer. For example 
there is no hope to locate critical points of the symbol starting only form 
the spectrum. 

The results we would like to present are perhaps of particular interest for 
^.-quantized Schrodinger operators on L 2 (R™): 

P h = -h 2 A + V, 

also called semi-classical Schrodinger operators. Here we will assume that 
the potentials V are smooth on R n and bounded from below. For this class 
of operators, the question is then to understand how certain fluctuations in 
the semi-classical spectrum can describe the shape of V. But our results will 
also apply to more general operators like /i-pseudo-differential, /i-admissible 
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operators (see definitions below) or Schrodinger operators on a Riemannian 
manifold M (simply replace the Laplace operator by a Laplace-Beltrami op- 
erator Am, see section 5). Most of these modifications are possible because 
our methods are micro-local and do not use global results on these opera- 
tors. 

Notations. Before entering into the details, let us give some definitions 
and recall some basic facts about the spectral theory for Ph- By a standard 
result, see [3], when V is bounded from below and with tempered growth, 
Ph has a self-adjoint realization on a dense subset of L 2 (M. n ). As usually, 
to this quantum operator Ph we can associate a classical counterpart with 
the Hamiltonian function p(x,^) = £ 2 + V(x), or total energy, on the phase 
space M. n x M n . In what follows, we note $f the flow of the Hamiltonian 
vector field: 

Hp = d^p.d x - d x p.d£. 
A classical energy surface is: 

= {(x, f ) G R n x R" / f + V(x) = E} 

and for a classical energy E GM we will say that: 

• E is regular when dp ^ everywhere on Eg. 

• E = E c is critical if dp = somewhere on T,e c - 

In the present article we are very precisely interested in a relation between 
the asymptotic properties as h — > + of eigenvalues Xj(h) of Ph : 



and the set of fixed point for the map &t (viewed as a map on M.t x T*M ?1 ). We 
recall in section 2, a sufficient condition to get discrete spectrum. It is well 
known, see [11] . that this semi-classical problem and the high-energy limit 
for the spectrum of the Laplace operator Am > on a compact Riemannian- 
manifold M are related. This relation can be viewed by quantizing 1/h = 



where Xj —> +oo is the increasing sequence of eigenvalues of Am. 
But in this article we will mainly consider the semi-classical problem and 
associated micro-local methods. 

Duality between the quantum and classical worlds. In geometry 
spectrum and periodic orbits can be related, in a very explicit way, by means 
of the Selberg (e.g. see [22]) and Duistermaat-Guillemin |14j trace formulae. 
It is in general in the most atypical situations, compact surfaces of constant 
negative curvatures or at the opposite completely integrable systems (e.g. 
like the free Laplacian on a flat torus) that the most explicit and exact results 
can be obtained. In quantum mechanics, the existence of such a relation is 
strongly suggested by the correspondence principle which asserts that, in 
the semiclassical regime h — > 0, many properties of Ph can be related to 
integral curves of &t and many invariant attached to the flow around these 
curves (see Eq.(pQ) below). For a general Hamiltonian H, not necessarily 
of the form kinetic energy plus potential, this correspondence principle is 



P h (fj(x,h) = Xj(h)<fj(x,h), Lfj G L 2 (M. n ) 
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also true under very reasonable assumptions on the symbol h(x, £) of H, a 
function on the phase space. 

In physics a more precise formulation of this principle appeared in the 
works of Balian&Bloch [2] and Gutzwiller |18j . The Gutzwiller formula is 
usually written as a trace formula for the resolvent of Ph at a given energy 
E: 

m V" 1 - Lvol (£g) i 1 V- a As-i 

(1) 2^x j (h)-E--^hyr + ih ^^ eh ' 

where in the r.h.s the sum concerns the closed orbits 7 inside T, E - Here 
Lvol(E£) is the Riemannian volume of T,e (defined w.r.t. the invariant 
Liouville measure), 5 7 = J £dx and Ay are respectively the classical action 
and the stability factor (including the Maslov phase) of the curve 7. Recall 
that the Liouville-volume LvoI(Se) satisfies the co-area formula: 

J Lvol(T lE )dE = J dxd£ = Vol^n^ 1 ([a, b]), 

[a,b] {(x,():a<p(x£)<b} 

where the measure on the r.h.s. is the Lebesgue measure of the pull-back. 

In mathematics and in physics, such a relation between spectrum and pe- 
riodic orbits provides a powerful tool of analysis and computation. See e.g. 
|25j concerning the asymptotic behavior of eigenvectors <Pj(x, h) and [19] for 
various applications in quantum chaos. See also [11] for a nice overview and 
applications in Riemannian geometry. 

Mathematical problems. For a Schrodinger operator on M. n , it is easy to 
check that two different type of divergence generally occur in Eq.(TjQ) : 

1) The sum over the spectrum is divergent when the resolvent is not a trace- 
class operator. In particular this is the case if V does not go fast enough to 
00 when |x| — > 00. When the sum appears to be convergent it can also have 
a divergent behavior when h — > 0. Worst, it can be that both sides of Eq. 
([T]) do not fit in the regime h — > 0. 

2) The sum over closed orbits is generally divergent. This is the case if |A-J 
does not decrease fast enough or if the number of periodic orbits of period 
smaller than T is exponentially growing with T. 

For example, if n = 1, using scaling and the asymptotic properties of the 
spectrum (here simply given by some Bohr-Sommerfeld quantization condi- 
tions, see [3]) it is easy to check that the trace of the resolvent of : 

d 2 

Q h = -h 2 -—r + \x\ a , h > 0, a > 0, 
dx z 

exists if and only if a > 2. Here, the harmonic oscillator, obtained for a = 2 
for which Xj(h) = h(2j + 1), is the limit case and the series diverges like the 
harmonic series. 
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1.2. Mathematical approach of the Gutzwiller formula. 

As seen above, the question to remove divergences has many important 
implication^] and we explain now a mathematical way to solve this problem 
via a smoothing of the so-called spectral density. To begin the discussion, 
simply assume that: 

For some E £ M, the spectrum of Ph is discrete, with finite 
multiplicities, in the interval [E — e, E + s], e > 0. 
A sufficient condition to obtain this property, uniformly w.r.t. E, is given 
in section 2. A well-posed problem is to study the asymptotic behavior of 
the spectral distributions: 

(2) T(E, h,<p)= <P( Xi{h) h ~ E ), as h -> 0. 

\Xj(h)-E\<e ' 

where ip is a test function chosen to remove the divergences. We can justify 
this terminology if we observe that the truncated spectral distribution: 

T E ,e{x) = 6 x 3 (h){x), (8xo, f) = f(x ), 

\\j(h)-E\<e 

acting on a function (p shifted by E and scaled w.r.t. h provides : 

,x-E' 



1{E,h,<p) = (T E>e (x),ip( h 

In reality this scaling w.r.t. h is very important and is used to get para- 
metrices involving the classical dynamics. Also it is not very hard to verify 
that when ip G <S(M) the size of the truncation (materialized here as e) is 
irrelevant on a scale of size 0(h°°) as long as e is strictly positive. I refer 
to section 3 for these points but I simply recall that 0(h°°) is the class of 
functions of h being in 0(h k ) for every k € N near h = 0. 

Statistical quantum mechanics. In general, apart in some very specific 
situations, it is not possible to compute explicitly the spectrum of Ph and a 
motivation to do semi-classical or high-energy estimates is to derive statistics 
about eigenvalues and their distribution. For example, in Eq.([2|) the formal 
choice of ip as the characteristic function of [— 77,77], < r/ < e, determines 
the number N(E, h) of bound states in [E — rjh, E + rjh]. 

This formal correspondence between T and the micro-local counting func- 
tion N has a mathematically rigorous formulation in term of Tauberian- 
theorems, see e.g. 0]. Under certain (generic) conditions^ on the symbol 
p of Ph, it can be proven that N(E,h) is proportional to /i 1_n times the 
Liouville- volume of the energy shell T,e- 

N ( E > h ) ~ fo Mn 1 LvQl (£g) aS h -»• 0+ - 

[2irn) n 1 



^Many important questions concerning the range of trace formulae (e.g., their validity 
beyond the Ehrenfest-time) are still open. We do not discuss these questions in this article. 
2 In particular the condition that E is non-critical for p, see below. 
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This is a micro-local formulation of the Weyl-law. A fortiori when n > 2 this 
implies that the finite sum defining T(E,h,<p) will involve a large number 
of eigenvalues in the regime h — > + . In general the formula for N(E,h) 
can be formally integrated to obtain a formula for counting eigenvalues in a 
compact interval: 

M([a,b],h) = #{j G N : Xj (h) G [a,b]} ~ VoW^M). 

Of course if the operator is bounded from below you can also use M(x, h) := 
M(] — oo, x], h). The question to estimate the remainder function for N is in 
general a relatively complicated problem and requires to use the properties 
of the underlying classical dynamics inside p _1 ([a, b]). Several other related 
problems, Riesz-moments or Lieb-Thirring inequalities, can be formulated 
in terms of N(E, h) and these problems are important in the 'stability of 
matter' problem. See [26] for an overview and references. Finally, I men- 
tion that certain Schrodinger operators with very singular critical sets (e.g. 
see [HE]) or non-confining potentials (e.g. see [30]) can lead to some very 
different kind of 'Weyl-asymptotics' for N or N . 

Relation with the classical dynamics. In reality the quantity T(E, h, ip) 
contains many interesting information (a priori more than the counting func- 
tions) since the asymptotic expansion of T(E,h,(p) when h — > + involves 
explicitly the classical dynamics on Eg and in particular the set of fixed 
point of the flow inside the energy surface: 

Fix£ = {(T,x,0 £lxE E : *t(x,£) = (x,0>- 

We recall that E is regular if Vp(x,£) / on and critical otherwise. 
Every critical point (#o,£o) ^ of p is a fixed point of our flow $ 4 since 
H p {xo,^o) = 0. In this situation we have R x {(xo,£o)} c Fix£ c . 

When E is not critical and the periodic orbits satisfy a condition of non- 
degeneracy, the asymptotics behavior of Eq.([2]) is well determined by the 
closed orbits of $j on Eg and the geometry of Eg. For the full treatment of 
this problem, and a complete formulation of the asymptotic expansion, we 
refer to [3 [27]. 

Removing divergences. We explain now shortly why the problem stated 
in Eq.([2]) leads to a mathematically rigorous version of the Gutzwiller for- 
mula. First, for each h > the sum is finite and a fortiori convergent. A 
convenient choice of tp also ensures that this quantity has an asymptotic 
expansion when h — > independently from the choice of e > up to correc- 
tions of order 0(h°°) as long as e stays strictly positive. Such a difference 
of size 0(h°°) plays no role since the discussion will be based on some finite 
order asymptotics w.r.t. h. 
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On the other side, only the periods of &t inside supp((/3), the support of 
the Fourier transform: 



contribute in the asymptotic expansion. This principle is useful since when 
supp((^>) is compact then finitely many closed orbits of Y^e contribute and 
the second divergence is solved. Hence if (p G Co°(R), the space of smooth 
functions with compact support, (p is in the Schwartz space Since ele- 

ments of 5(IR) are smooth with exponential decay at infinity, no divergence 
occurs and the size of e is irrelevant, up to an 0(/i°°)-error. 

Finally, in Eq.([2]) the scaling w.r.t. h is important. With this choice and 
via Fourier transform considerations, we can use the semi-classical propaga- 
tor Uh(t) = exp(itPh/h), solution of the Schrodinger equation: 



to obtain a precise control w.r.t. h. Roughly, Uh(t) can be expanded w.r.t. 
h via a so-called WKB approximation. This expansion also provides the 
explicit relation with the classical dynamics. The precise technical justifica- 
tions are given in section 3. 

1.3. Critical values and contributions of equilibria. In the previous 
section we heuristically outlined a relation valid when the semi-classical 
parameter tends to 0: 



Meaning that the asymptotic behavior of the left hand side can be expressed 
in terms of distributions generated by fixed point of the flow. 

In the r.h.s any point (x, £) of a periodic orbit 7 appears only at times 
kTj, £ Z*, where T§ is the primitive period of 7 orbit. Also any point 
of the energy surface contributes for t = since the flow is the identity at 
t = 0. But an equilibrium (xq,£q) satisfies <3?j(xo,£o) = (^o>£o) for all t. 
Hence when E is no more a regular value the nature of the set of fixed point 
changes and some new contributions appear in the asymptotic expansion. 
These new contributions can be qualified of new wave invariants (see below) 
and are extremely important for the inverse spectral problem. 

When E = E c is a critical value of the principal symbol p, the asymptotic 
behavior of Eq.([2]) is more complicated and is closely related to the geometry 
of the flow inside £_e c . The presence of classical equilibria inside Y>e c and 
the stability of the flow near the critical set affects strongly the nature of 
the asymptotic expansions. For a non-degenerate critical point, i.e. when 
d 2 p(xo, £0) is an invertible matrix when dp(xo, £0) = 0, the reader can consult 
[4j. The problem is treated there for quite general operators, also including 
the case of a manifold of critical points, but for supp(</?) small around the 
origin. For Schrodinger operators on lR n and supp(</?) compact but arbitrary, 
the results of [4] are improved in [24]. 



<p(t) 




ihd t u h (t) = p h u h (t) 



lim 1{E,h,<p) ^ Fix £ = {(t,x,£) £lxS E / $ t (x,0 = {x,§}. 
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Two important problems occur in presence of critical points. First, at 
every point where dp = the surface T,e c and the metric of T<e c are not 
smooth. Next, the determination of the asymptotic expansion w.r.t. h can 
be very difficult. The point is that T(E, h, (p) can be expressed in terms of 
oscillatory integrals: 

1(h) = J a(t,x,Oe^ {t ^°dtdxdt h^0 + . 

MxT*R n 

The oscillating coefficient h^ 1 is precisely imposed by the scaling w.r.t. h 
in Eq.(|2]) an d plays an important role since 1(h) oscillates fast in the semi- 
classical regime. Via the WKB approximation, the phase if) is related to the 
flow so that the asymptotic behavior of 1(h) is determined by the closed 
orbits. The technical problem is that, in presence of an equilibrium, tp has 
some degenerate critical points. The stationary phase method cannot be 
applied and the asymptotic expansion of 1(h) is radically different : e.g. 
some terms h a , a £ Q and powers of log(/i) generally appear in this setting 
(see below). Also the nature of these new terms can be very different since 
for example they can be associated to some distributions acting on (p with 
a continuous support (e.g. the full set of real numbers or a half-line). 

Wave- invariants. A classical approach (used in section 3) is to study the 
asymptotic behavior, as h — > + , of the localized trace: 

n(E,h,t) = Tr (e(P h )e~^ Ph - E ^ ,0 € . 

I will follow now the terminology used in [23] • Under certain assumptions 
(see section 3 and 4), and for E regular it is well known that 0, admits an 
asymptotic expansion of the form: 

oo 

j=-n 

The coefficients a,j(E,t) are some tempered distributions on the line Mt 
and are called wave invariants of Ph. These distributions have a different 
expression when E varies and many of them are continuous functions of E 
as long as we do not cross critical levels of the energy function p. 

When E — >■ E c (E c stands for critical levels) one can observe a discon- 
tinuity in the asymptotic expansion. Also some new coefficients generally 
appear, since we have the asymptotics: 

n— 1 oo 

(3) n(E c , h, t) ~ ^2 ^2 a jA E c t)h* log(h) k , as h -»■ + , 

k=0j=-n 

for some p G N*, see [H El [H [91 [23] for details and examples. These new 
coefficients are called new wave invariants and the top order coefficient of the 
expansion in Eq.( [3"]) contains many information on the shape of the symbol. 
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1.4. Results and strategy. Our first objective is to relate some variations 
in the discrete spectrum of Ph with the presence of fixed points for the 
classical system: this principle detects the presence of new wave invariants 
and a fortiori critical energy levels. Secondly, we establish that the precise 
knowledge of such a spectral fluctuation can describe the singularity of the 
potential (or of the symbol for general operators). In theory, such a determi- 
nation is possible since the contributions of equilibriums are highly sensitive 
to the local shape of V and are extremely persistent when the test function 
<p varies. 

For Schrodinger operators, we will consider the case of a potential V with 
finitely many critical points Xg attached to local homogeneous extremum of 
V. An immediate consequence is that p admits, locally, a unique critical 
point (4,0) on the surface = {(x,£) G R 2n / £ 2 + V(x) = V(x J Q )}. A 
typical example is a polynomial double well in dimension 1 where 3 critical 
points occur at the 2 minima and at the maximum of V (see figure 1). 




Figure 1. Non-symmetric and symmetric double well. 



In fact using certain generalizations of stationary phase methods, neces- 
sary if the phase has some degenerate critical points one can derive a very 
precise relation between the spectrum and the set of fixed points of the flow 
also including the new wave invariants. Once this relation is established in 
the form of asymptotic expansion, the main results follow since: 

• Equilibriums have a continuous contribution w.r.t. the time t. 

• Shrinking supp((^) erases all other contributions: Weyl-terms and 
periodic orbits. 

• The remaining contribution, given by new wave invariants, displays 
some nice information on V. 

The first assertion simply means that a fixed point generalljH contributes 
to the asymptotic expansion of T(E J c ,h,ip) in the form h a log(h)P (T at p,tp) 
where T a ^ is a distribution such that supp(T Q)i g) = R, [a, oo] or [— oo, a]. 
Contrary to standard periodic orbits whom contributions are supported in 
the set of periods, such a term supported on the line cannot be erased just 
by shrinking the support of <p. 



Some discrete contributions can also sometimes occur as pointed out in [4] or [6]. But, 
to attain our objectives, we can avoid to include them in the spectral estimates. 
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For example, if supp(</?) contains no period of the flow our analysis follows 
if we view T(E, h, (p) as a function of E : 

• The order w.r.t h of T(E, h, p) changes when E — >■ E c (Prop. [5j). 

• This discontinuity at E J C describes the shape of V. 

This indicates the presence of an equilibrium for &t, a fortiori of a critical 
point for V. 

Remark 1. For a degenerate singularity the information is more difficult 
to interpret compared to a non-degenerate singularity (see section 2). 

2. Hypotheses and main result. 

Let p(x, £) = £ 2 + V(x), where the potential V is real valued and smooth 
on M. n . To this Hamiltonian is attached the operator Ph = —h 2 A + V(x) 
and by a classical result, see [3], Ph is essentially self-adjoint starting from a 
dense domain of L 2 (R n ) when V is bounded from below and with tempered 
growth. 

Remark 2. In this section, we are here mainly interested in the case of 
Schrddinger operators but generalizations to an /i-admissible operator (e.g. 
in the sense of |29j ) are given in section 5. 

First, to obtain a well defined spectral problem, we use: 
(Hi) V £ C 00 (W n ). There exists C £ M such that liminf V > C. 

oo 

Note that (Hi) is always satisfied if V goes to infinity at infinity. Now, con- 
sider an energy interval J = [E\,E<2\ with E2 < liminf V. In the following 
we note : 

(4) J{e) = [Ei-e,E 2 + e\. 

For e < £0 the set p~ l (J(e)) is compact. By Theorem 3.13 of |29j the 
spectrum a(Ph) H J(e) is discrete and consists in a sequence: 

X x {h) < X 2 (h) < ... < Xj(h), 

of eigenvalues of finite multiplicities, if e and h are small enough. In general 
such a condition that the pullback of J or J(e) by p is compact is sufficient to 
obtain a discrete spectrum. This is not necessary as shows the non-confining 
potential x 2 y 2 on M 2 , see [30J where precise spectral estimates are given for 
such potentials. 

The central object of study is the spectral distribution: 

,Xj(h)-E 



(5) ?(E,h,<p)= V{ 



h 

Xj{h)&J{e) 

and, more precisely, the asymptotic information contained in this object as 
h — > + . To avoid any problem of convergence we impose the condition: 

(W.2) We have <p G Co°(M) with a sufficiently small support near the origin. 
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Remark 3. (H.2) is used to erase contributions of non-trivial closed orbits 
and can be relaxed to (p G Co°(K) with a weaker result. A more precise de- 
scription of supp(i^) is given in Lemma [TBI For a non-degenerate minimum, 
it is more comfortable to assume that supp((^) contains no period of d$t(^o). 
Some singularities, not directly relevant here, are generated by these periods 
and we refer to [U [24] for a detailed study of these contributions. 

To simplify notations we write z = (x,£;) G M? n and Y*e = P ^({E}) and 
we use the subscript E c to distinguish out critical values of p. Of course one 
can also work with T*M n ~ M n x M n . In J there is finitely many critical 
values E\ , E l c and in J) finitely many fixed points z$, Zq 1 , m > I. 
We impose now the type of singularity: 

(W3) On each S^j the symbol p has isolated critical points Zq = (xq,0). 




These critical points can be degenerate but are associated to a local extremum 



where V^k, homogeneous of degree 2k, is definite positive or negative. 

Remark 4. For non-degenerate singularities we can apply the results of 
[H El [23] and the extremum condition is not really necessary. We will recall 
and use these results in the next section. But for a degenerate critical 
point of V the extremum condition is required since, to our knowledge, the 
contribution of such a singularity to the wave expansion is unknown. 

The next assumption, erases the mean values, i.e. the heat-like invariants 
or so-called Weyl-terms, in the trace formula: 

{U4) <p is flat at 0, i.e. <p^{0) = 0, Vj G N. 

We can weaken condition (^4) to <^')(0) = 0, Vj < ko, where ko G N* 
depends only on the degree of the singularities of V (see section 4) without 
essential change. Such a function <p exists and is easy to construct. Pick <j> G 
C$°(M), supp(0) C [-M,M], then <p{t) = t 2 i°(t>{t) satisfies our hypotheses. 
In this case, we can chose the function (j) even so that ip is real. 

Finally, to relax a bit (H.2) we need a control on the contribution of closed 
orbits. To do so, we impose the classical condition : 

(H5) All periodic trajectories of the flow are non- degenerate. 

Non-degenerate closed orbits are those whose Poincare map does not admit 
1 as eigenvalue and are isolated. The non-degeneracy condition on orbits is 
not a central argument in this work and is only used to control the order 
w.r.t. h of the contribution of closed orbits (i.e. the order of the usual 
wave-invariant). One could also impose a condition of 'clean- flow' to con- 
sider families/submanifolds or bunches of closed orbits of positive dimension. 
These conditions on the classical dynamics can be simply discarded as soon 
as we have a strictly positive lower-bound on periods of closed-orbits like in 
Lemma [16J 




ofV: 
(6) 
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Detecting critical levels. The first result shows how to detect critical 
energy levels by revealing a singularity in the spectral estimates: 

Proposition 5 (Spectral variation). 

Assume that conditions (Hi) to (Ha) are satisfied. As h tends to + , we 
have: 

T(F u n sj ( h °°) ^fEG[E 1 ,E 2 }\{E 1 c ,...,E l c }, 

wtw-X oum i/E=Ei, j€{w>, 

where each fj(h) has a finite order w.r.t. h. 

For a non-degenerate critical point the coefficients fj can be determined 
explicitly It is in general possible to predict a full asymptotic expansion 
but to obtain an invariant formulation of all the coefficients can be difficult 
since for degenerate critical points the method you have to use is more 
complicated than the usual stationary-phase formula. 

For example, if E = E\ and the surface T, E j carries a single minimum of 
degree 2k, k > 1, using the main result of [8] we obtain: 

(7) f j (h) = D(n,k,<p)h^- n . 

In the same situation, but for a local maximum of V, using the results of 
[5], we can obtain a logarithm of h: 

(8) fj(h) = D(n,k,ip)h% + %- n log{h) j , j = or 1. 

In Eq.(J7]) and Eq.(|8]) the coefficient D is a tempered distributions acting ip 
characteristic from the nature of the critical point. In fact if the critical 
surface carries more than one critical point then fj is the sum of their 
respective contributions. Note that for n = 1 and k > 1 the singular term 
has negative order w.r.t. h. A detailed formulation of the coefficients fj(h) 
is given in Propositions I18|19l and [20l 

Inverse spectral results. An interesting property is that in the singularity 
of T(E, h, ip) when E -4 E c the order w.r.t. h but also the constants of the 
top order-coefficients, describes partially the shape of V: 

Theorem 6 (Inverse result for Morse-critical points). 

Assume that Y<e c carries exactly one critical point associated to a non- 
degenerate critical point xq of V . Then the discontinuity of T(s,h,(p) at 
s = E c determines the spectrum of d 2 V(xo). 

This result follows from the special form of the Duistermaat-Guillemin- 
Uribe density at a non-degenerate critical point. Observe that in particular 
we retrieve the Morse index of V at xq (number of positive eigenvalues 
minus number of negative eigenvalues of d 2 V(xo)). Once more, because of 
the invariance under coordinates permutations, or under a rotation of the 
potential around the critical point xq, it is in general not possible to retrieve 
the quadratic form d 2 V(xo) in a given system of coordinates x. 

For a degenerate homogeneous singularity we have also a nice result: 
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Theorem 7 (Inverse result. Degenerate critical points of V.). 

Assume that there is exactly one critical point (xq,0) on the singular energy 
surface Tie c - Assume that xq is attached to an homogeneous maximum or 
minimum of the potential of degree 2k. Then the discontinuity ofT(E, h, (p) 
at E = E c determines: 

• The degree 2k of the critical point ofV. 

• The spherical mean-value of the germ ofV in xq: 



Observe that A(V) itself is invariant under rotation and translation of 
V . Both results of Theorems [6] and [7] are limited in presence of multiple 
equilibriums on the same surface since the sum of contributions of each crit- 
ical point could lead to a compensation or to several contributions of exactly 
same nature and order. In general if V is not a Morse function, or if a surface 
of energy T,e c carries more than one critical point, eigenfunction-estimates 
seems to be required to have a well-posed inverse problem. Finally, at the 
end of the article we will show up some interesting invariants for pseudo- 
differential operators with homogeneous singularities attached to extremum 
of the symbol. 

Remark 8. I would like to emphasize that a maximum is more difficult to 
detect contrary to a local minimum which is an isolated point of the energy 
surface (locally £e c is just a point). A similar result holds for an operator 
of the form T(£) + V(x) (kinetic plus potential energy) where T is convex 
near the origin and V has a local minimum at xq. Moreover a maximum, 
attached to an unstable critical point of the flow, is much more complicated 
to treat with semi-classical methods. See [9] for a detailed study. 

Longer range estimates. In (H2) the condition that supp(c/3) is small 
implies a very accurate spectral estimate (e.g. by a Paley- Wiener estimates 
for the decay of <p). It is possible to relax this assumption but the result is 
a bit weaker: 

Corollary 9. Assume that conditions ('Hi), (%), (H4) and (H5) are sat- 
isfied and that (p E Cg°(M), then we obtain: 



For critical values of p, estimates are the same as in Proposition^ 

The justification, see section 4, is that in this case the asymptotics is given 
by a finite sum over periodic orbits of energy s. This result is weak if the 
singularity of V is non-degenerate since the equilibrium has a contribution 
of degree w.r.t. h, see Propositions I18|19|20l or section 3 of [4]. 

In theory there is always a variation when E —¥ E c but of course this 
effect can be harder to detect if there is no change in the order w.r.t. h. In 



A{V) 




T{E,h,<p) =0(1) ifEe[E l ,E 2 ]\{E 1 c ,...,E l c }. 
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that situation there is only a discontinuity in the top-order coefficient w.r.t. 
h so that the result can be qualified of 'weaker'. 

3. Oscillatory representation. 

The construction below is more or less classical and will be sketchy. The 
only change with the usual construction, around a single energy level, is that 
we use a more global localization around J = [E\,E<i[. Strictly speaking, 
with (Hi), we could also consider ] — 00,^2] since there is no eigenvalue 
below a fixed energy level Eq given by the minimum of the quadratic form 
attached to P^. Let be <p G <S(R) with tp G Co°(R), we recall that: 

T(E,h,<p)= Yl <P( Xj{h l~ E )> J(e) = [Ei-e,E 2 + e], 

with p _1 (J(e)) compact in T*R n . For e > small enough, we localize 
around J with a cut-off G G Cq°(]Ei — e,E2 + e[), such that = 1 on J 
and < < 1 on R. We accordingly split-up our spectral distribution as: 

T(E, h, ip) = Ti(E, h, p) + T 2 (E, h, 99), 



with 



Ti(E,h,<p)= (l-G)(*j(h)M Xj{h) h E ), 

Xj(h)eJ(e) 

r 2 (E,h, V )= y, ®(Wh)M Xj{h) ~ E y 

Since (p G <S(R) a classical estimate, see e.g. Lemma 1 of [7], is: 

(9) Ti(E, h, ip) = 0(h°°), as h -»• + . 
By inversion of the Fourier transform we have: 

Q(P h ) (p (^l) = LJ e ^p(t)exp(- l -l-P h )e(P h )dt. 

R 

The trace of the left hand-side is T 2 (E, h, tp) and Eq.@ provides : 

(10) T(E, h, <p) = ^Tr j <p(t)exp(~P h )e(P h )dt + 0(h°°). 

R 

Eq. ()10p is very close to the classical Poisson summation formula on S , see 
|31j for a discussion and an interpretation, since the r.h.s. is expressed 
below in term of the classical dynamics. This asymptotic relation justifies 
the terminology of trace formula. 

Moreover, the formulation in Eq. (|10p shows that the scaling w.r.t. h, 
imposed in the initial definition of T(E,h,(p), is the best one since we will 
solve now the semi-classical propagator homogeneously w.r.t. h. Let Uh(t) = 
exp(— t^-P/i) be the quantum propagator. We approximate Uh(t)Q(Ph) by a 
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Fourier integral operator (FIO) depending on h. Let A be the Lagrangian 
manifold associated to the flow of p: 

A = {(t, r, x, £, y, V ) e T*R x T*R n x T*R n : r = p(x, £), (x, £) = * t (y, ??)}, 

and I(R 2n+ , A) the class of oscillatory integrals based on R 2n+1 and whose 
Lagrangian manifold is A. The next result is a semi-classical version of a 
well known result on the propagator, see e.g. Duistermaat |13| . 

Theorem 10. The operator Uh(t)Q(Pf l ) is an h-FIO associated to A. For 
each N 6 N there exists U^^t) with integral kernel in Hormander's class 

I(R 2n+1 ,A) and R^\t) bounded, with a L 2 -norm uniformly bounded for 
< h < 1 and t in a compact subset of R, such that: 

U h {t)Q{P h ) = U^{t) + h N R[ N \t). 

This result provides the existence of an asymptotic expansion in power of 
h with a remainder that can be controlled since supp((/3) is a compact. After 

perhaps a reduction of e, this remainder (t) is estimated via: 

Corollary 11. Let Q 1 £ C^°(R), with Q 1 = 1 on supp(0) and supp(Gi) C 
}E 1 - 2e, E 2 + 2e[, then \/N € N: 

Tr(e(P fe M^^)) = ^Tr / m^ tB U^(t)Q 1 (I) i )dt + 0(h N ^). 

R 

For a proof of this result, based on the cyclicity of the trace and a priori 
estimates on the spectral projectors (see we refer to |7]. For the 

particular case of a Schrodinger operator the BKW ansatz shows that the 

integral kernel of Uq^ (t) can be recursively constructed as: 

M" 

b[ N) = b + hb 1 + ... + h N b N , 

where S satisfies the Hamilton-Jacobi equation: 

p{x,d x S{t,x,Q) + dtS{t,x,O=0, 
with initial condition S(0, x,£) = (x,£). In particular we obtain that: 

{(t,d t S(t,x,r}),x,d x S(t,x,rf),d ri S(t,x,ri),—7])} C A, 
and that the function S is a generating function of the flow, i.e.: 

(11) Md v S(t, X, 7?), 7?) = (x, d x S(t, X, 7?)). 

We insert this approximation in Eq. (|10p . we set x = y and we integrate 
w.r.t. x. Modulo an error 0(h N ~ n ), we obtain that T(E,h,ip) equals: 

(12) — ^ / eW^-<'^a^(t,x,mt)dt^, 

MxT*M" 
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where a h {t,x,rj) = b h (t,x,x,r]). 

Remark 12. By Theorem 3.11 & Remark 3.14 of [29], G(P/ l ) is /i-admissible. 
Moreover, the symbol is compactly supported in p~ l {[E\ — e, E2 + s]). This 
global result w.r.t. E £ [E\,E2\ allows to consider below only oscillatory 
integrals with compact support for the evaluation of the spectral distribu- 
tions. 

Micro-localization of the trace. 

If ip £ C^°(T*M n ), we recall that ip w (x, hD x ) is the linear operator obtained 
by Weyl-quantization of i/j and semi-classical quantization of £. This means: 

R 2n 

Mainly, the contribution of an equilibrium zq £ Y<e c can be reached via: 

(13) T Z0 (E C , h, ip) = i-Tr J e^m^ix, hD x )exp(~tP h )e(P h )dt, 

R 

where ip £ C^(T*M n ) is equal to 1 near zq. Since the trace is a cyclic 
operation, we can use the symbolic-calculus and to insert an L 2 -bounded 
observable (here a cut-off in the phase space) is in general not very expen- 
sive in the FIO construction. This construction with smooth cut-off allows 
to work with pseudo-differential partition of unity. This approach is use- 
ful to obtain a weak generalization of our results in presence of multiple 
equilibriums. 

We recall some basic results on the symbolic calculus with FIO. Horman- 
der's class of distributions with Lagrangian manifold A over M n is noted 
I(R n ,A). If (x ,£q) £ A and (p(x,6) £ C°°(R n x R N ) parameterizes A in 
a sufficiently small neighborhood U of (xo,£o), then for each Uh £ I(R n ,A) 
and x ^ Cg°(T*IR n ), supp(x) C U, there exists a sequence of amplitudes 
Cj(x, 9) £ C£° (R n x R N ) such that for all L £ N: 

X w (x,hD x )u h = hn( Cj e^) + 0{h L ). 

-d<j<L 

Hence, for each N £ N* and modulo an error 0(h N ~ d ), the localized trace 
T Z0 (E C , h, if) of Eq. (fT3"|) can be written as: 

(14) (2nh)- d j eS^'O-^-O+^a^^t.x.O^Wdtdsde. 

RxR 2 ™ 

To obtain the exact power — d of h, we apply results of Duistermaat [13] on 
the order of FIO. Since: 

• /i-pseudo-differential operators ip w (x, hD x ) are of order w.r.t. 1/h. 

• the order of U h (t)@(P h ) is -\. 
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if we identify the operator with it's distributional-kernel we have: 



Multiplying by (p{t)eh tbjc and passing to the trace we find Eq. (|14p with d = n 
and we write again d^\t,x,7j) for the diagonal evaluation d^\t, x, x,rj). 



In particular: 

(15) af£\t,x,x,rj) = ip(x,r})a (t,x,x,r)), 

is independent of h and is compactly supported w.r.t. (x,r]) since V £ 

4. Proof of the main result. 

Let be E c any critical value in [E\, E2} and zq an equilibrium of T,e c - We 
choose a function ip G C^°{T*W n ), with tp = 1 near zq, hence: 

T 2 (E c ,h,<p) = ±-Tr J e itJ tm^ W (x,hD x )eM-{tP h )e(P h )dt 

R 

+ ^Tr y e itJ t<f(t)(l-i; w (x,hD x ))e W (-^tP h )9(P h )dt. 

R 

If there is no other singularity on T,e c with (^5) the asymptotic expansion 
of the second term is given by the semi-classical trace formula on a regular 
level. For finitely many critical point on E£ c , we can repeat the procedure. 
The first term is micro-local and precisely generates the singularity in The- 
orem We note Q the discrete set of critical points Zq in p _1 (J). 

Classical dynamics near the equilibrium. 

A generic critical points of the phase function of Eq. (fT2j) satisfies the rela- 
tions: 

The right hand side coincide with the set Fix^ of fixed points of <&t defined 
in section 1. This set generally consist of closed trajectories of the flow 
inside Eg, the energy surface (for t = the flow is the identity) and finally 
equilibria/critical-points. By the non-stationary phase lemma, outside of 
the critical set Fix^ the contribution is of order 0(h°°). 

To apply the stationary-phase methods it is important to study the nature 
of the phase function along the critical-set Fbcg. For the regularity of the 
Hessian of the phase function, the next lemma is particulary useful. We will 
often denote points 6 T*R™ of the phase space by a single letter z. 
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Lemma 13. Let us define ty(t,x,£) = S(t,x,£) — +tE c , then if zq is 
critical point of ^ , we have the equivalence: 



In other words, degenerate directions of the phase correspond to fixed points 
of the linearized flow at zq . 

The proof is standard and can for example be found in [Ml E]. We recall 
that the linearized flow d§t is the differential of the flow w.r.t. initial 
conditions z = (x, £). If we use Lemma [T3| we obtain for our phase function 

Corollary 14. A critical point (T, xo>£o) °f ^(ti x i£) * s degenerate with re- 
spect to (x,£) if and only ifT is a period of the linearized flow d x £&t(%0) Co)- 

The next result is also well known, see, e.g., [lj, from classical mechanics 
and differential geometry: 

Lemma 15. If d x p(x ,£ ) = %>(xo,£o) = then d X) £$ t (? , Co) is the 
Hamiltonian flow of the quadratic form ^d 2 p(xo, £o) on T XQ ^ (T*M. n ). 

Hence, when zq is a critical point of p, the linear map w i— > d$t{ z o) w 
can be interpreted as the Hamiltonian flow of u i— > Q(u) = (^d 2 p(zo)u,u) . 
Observe that for a Schrodinger operator a critical point is always of the form 
zq = (xq,0) with dV(xo) = and the quadratic form is: 



A fortiori the map d§t(zo) is governed by a quadratic Schrodinger operator. 
Degenerate critical point. If the kernel of d 2 V(xo) is not trivial and 
contains a non-zero vector v the linearized-flow in the direction of v is the 
flow of the free Laplace operator and we have d& t (zo)(v,0) = (v,0) for all 
t. Hence it is never possible to apply the Morse-lemma, and a fortiori the 
stationary phase method (and this for any time t). But if we restrict the 
study to some homogeneous singularity for V it is possible to find a local 
diffeomorphism changing the phase into a polynomial function R (i.e. a 
local normal- form for For these normal forms it is possible to generalize 
the stationary-phase method: 



We will not review all these geometric and analytic results individually. We 
will refer to [U [TJ [U (H [10] for the reduction of the phase function and 
the asymptotic of the resulting oscillatory integrals. I mention that it is 
complicated to express all distributional coefficients Cj^kip) invariantly since 
a depends on the local-diffeomorphism x transforming ^ into R. The top- 
order coefficients are relatively easy to express, the next one start to depend 
on the derivatives of x- To express all Cj t k(a) in terms of V would improve 
the inverse spectral result to higher derivatives of V at the critical point. 



d 2 z V(z )6z = 0« d x $ t {zo)5z = Sz, Vfe e T Z() (T*R n ). 
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Non-degenerate critical point of V. When (PV(xq) is invertible the 
have a classical result. It is easy to check that the sign of the quadratic map 
Q determines the 'stable' and 'unstable' directions after diagonalizing d 2 V 
via an orthogonal linear transformation. 

Combining the results of Lemma [15] and Corollary 1141 to achieve our goal 
it will be sufficient to stay below the smallest positive period of the linearized 
flow. Working in some suitable local coordinates near zq we can assume that 
the quadratic form attached to the potential is of the form: 

r n 

V(x) = E c -^2ajXj + ^2 ajxj + 0(\\x\\ 3 ), ay > 0. 

j=l j=r+l 

This can be achieved by a translation, a change of linear coordinates (via an 
orthogonal matrix) and an eventual permutation of coordinates. Observe 
that all these linear transformations, in particular the action of the orthog- 
onal matrix, leave the Laplace operator invariant. Of course the variable 
attached to indices j = 1, ... , r will generate hyperbolic functions (and hence 
never degenerate in the sense of Corollary I14h . 

A lower bound on primitive periods. 

The next result provides a global information on the smallest positive prim- 
itive periods of the classical flow. This lemma will be used to extract the 
new wave invariants appearing at a critical energy level. 

Lemma 16. There exists T > 0, depending only on V and J = [Ei,E 2 ], 
such that $ t (z) / z for all z G p -1 ( J)\0 and all te]-T, 0[U]0, T[. 

Proof. If Hp is our hamiltonian vector field and z = we have: 

\\H p ( Zl ) - H p {z 2 )\\ 2 = 4||£i - 6|| 2 + ||V^(xi) - V x V(x 2 )\\ 2 . 

When z\ and z 2 are in the compact p~ 1 (J) there exists b > such that: 

\\V x V{ Xl ) - V x V(x 2 )|| < 6||xi - x 2 \\. 

Hence, there exists a > such that: 

H-Hp(zi) - H p (z 2 )\\ < a\\zi - z 2 \\, Mzi,z 2 £ p _1 (J). 

The main result of |32] shows that any periodic orbit inside p~ l {J) has a 
period r > 2ir/a > 0. The lemma follows with T := T(V, J) = lixja. ■ 

Remark 17. The result of [32] is optimal (for the harmonic oscillator the 
previous inequality becomes an equality). Note that T is decreasing if one 
increase the size of J. Lemma [16] provides a total control on the r.h.s. of 
the trace formula. If ip £ C^°(]0,T[), the only contribution arises from the 
set {(t, Zq), t £ supp(y?)}, i.e. from the new wave invariants. 

Now, we restrict our attention to the singular contribution generated by 
one critical point. We check now the condition of non-degeneracy of our 
phase-function. As it was explained in section 2, for a non-degenerate ex- 
tremum a minor technical problem could occur. If xq is a maximum of the 
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potential d$t( z o) nas no non-zero period which ends immediately the discus- 
sion. If xq is a minimum d$t(zo) is elliptic with primitive periods (T\, ..,T n ) 
generated by the eigenvalues of d 2 V(xo). But the constant b of Lemma [TBI 
is certainly bigger than the spectral radius of d 2 V(xo) and hence we have 
the inequality T < rnin{Ti, .., T n }. Following the approach of [H[B1[24|, if 
supp((^) c] — T,T[ the associated contribution is smooth on supp((^)\{0}. 

For a degenerate critical point zq as in (^3) a surprising result, established 
in [HI [9], is that the only singularity, for the first new wave invariant, is 
located at t = 0. Hence the condition <p £ Cg°(]0,T[) or <p £ C^(] - T,0[) 
is sufficient to determine the new- wave invariants generated by Zq. 

4.1. Non-degenerate critical points. With the previous considerations, 
we consider z$ E T*e c a non-degenerate critical point of p, a fortiori isolated, 
and we assume that supp(<£) c] — T, 0[U]0, T[. In this setting we know that 
the phase function ^ (introduced in Lemma [T3|) has a non degenerate critical 
point in zq for all t E supp(<£). 

In our setting, up to a change of local coordinates, we can assume that 
p(x,£) = p2(x,£) + 0(\\x\\ 3 (near the origin zq = (0,0)), where the quadratic 
form is: 



(16) p 2 (x,0 = - ajx*)) + ( Y, (C| + 

j=l j=r+l 

The flow of p 2 , viewed as an element of End(r (r*R n )) ~ End(R 2n ), is 



exp(tH P2 )(x,0 =A(t) 



A(t) 



( a{t) e(t) 

b(t) 

c(t) a(t) 

\ d(t) 



x 

\ 

/(*) 


b(t) J 

where (x,£) = [x' , x" , £"), x',£' E W, x",£" E M n - r , and 

a(t) = diag( ch(aji)), 
b{t) = diag( cos(ayi)), 
c(t) = diag( ajsh(ajt)), 
d(t) = diag( — ctj s'm(ajt)), 
e(t) = diag(^-sh(a i t)), 

f(t) = diag( J-sin(a^)), 

and the symbol "diag" means diagonal matrix. Observe that for a(t), c(t), e(t) 
the index j varies in {l,r}, for the others it varies in {r + l,n} and that 
detA(i) = 1 for all t. 

Local reduction of the phase function of our FIO. Since S is solution 
of the Hamilton-Jacobi equation we have (locally) the relation: 

*t(dtS(t,x,S),Z) = (x,d x S(t,x,0). 
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But since $t( z o) = z o> writing a Taylor expansion for: 

S(t,x,0 = S 2 (t,x,0 + O(\\(x,0\\\ 

and &t (always near the critical point zq = 0) we get that S 2 satisfies the 
relation: 

dM0)(d !; S 2l + O(\\(x,0\\ 2 ) = (x,d x S 2 ) + O(\\(x,0\\ 2 ), 

where the linear map cM?j(0) can be determined as above. This linear system 
is regular exactly when det d(t) ^ 0. As a consequence, we pick: 

2ir 

l = mi(—,je{r + l,...,n}), 

ay 

and an interval L c]0, l[. For each t € L, it follows that we can find a 
(time-dependant) change of coordinates (t, x,£) h-> (i,x), well defined in a 
small neighborhood of supp(c^) x {^o} such that: 

*(t,z)~Q(x), 

Q(x) = (iXl - X2) + • • • + (Xlr-l - Xlr) + (Xlr+l + X2r+ 2 ) + ' ' ' + (xL-1 + xL))- 

Applying the stationary phase method for the x variables, which is legal 
since the remaining integration w.r.t. t is of compact support, we get: 

T 20 (£ c ,^,/i)~eX^(Q) J 0(t)d fJ i t (z o ) + O(h). 

M 

Here sgn(Q) = (n + n — r) — r = 2(n — r), next by checking: 

• The value of x(t, z o), 

• The value of the Jacobian at the point zq, 

it comes out, % being given by x(i,z) = (d$t( z o) — Id))-2 + C(||z|| 2 ) at the 
first order, and as long as det(d&t(zo) — Id) / 0, that: 

(17) dutfa) = I det(d* t («b) - Id)|-5 

Observe that this coefficient, the Duistermaat-Guillemin-Uribe density, is 
indeed a smooth function as long as we stay away from any period of the 
linearized flow at Zq. An explicit computation using A(t) in our coordinates, 
done for example in [23], shows that the density is given by: 

dfh(zo) = — — • 

I J] smh(aj(z )t) J] sm(aj(z Q )t)\ 

j=l j=r+l 

On this formula we see that we have r-negative and n—r positive eigenvalues 
at the critical point xq. The desired result follows since d/j,t(zo) determines: 

• The signature of the Hessian of V at xq. 

• Eigenvalues ctj(xo). 
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The second affirmation follows by Taylor expanding the density and evalu- 
ating it at different points. 

Comments. 

Such a density was first introduced by Duistermaat and Guillemin [14] . But 
it seems that the exploitation of this term to describe contributions of crit- 
ical points goes back to Guillemin and Uribe |16j . This kind of density can 
be extended to Morse-Bott singularities for V (see [U El])- For a strict 
minimum of the potential xq the shape of cZ/it (xo,0) is fundamental to get 
better inverse spectral results. In particular one has to check the presence of 
resonant coefficients. See [12J and |17| for improved results near a minimum. 

These densities can be continued as meromorphic-distributions with sin- 
gular support: 

• At the origin, such a study is done in [3]. 

• With singular at a period T of d$t(zo), this is done in |24j . 

For an operator which is not a Schrodinger operator some new terms can 
generally appear at a period of d$>t(zo) (see [6]). All these facts strongly 
suggest that the 'pike singularity' of dfit(zo) as t — > T near a period T of 
d&t( z o) should also describe the symbol. For example, for a Schrodinger 
operator a double eigenvalue w > of d 2 V(xo) will generates a singularity 
of double magnitude in d/j,t(zo) at the point T = ^L. 

Example. An important toy model is the case of an n-dimensional har- 
monic oscillator: 



This model is one of the few Hamiltonians that can be explicitly solved. 
Then, for t smalfl and t ^ ^^7' the kernel of the propagator Uh{t) = 
exp(^P/j) can be explicitly computed as: 



This clearly shows that the small h behavior of Uh(t) determines the eigenval- 
ues Wk- This important example can be perturbed to use Dyson-expansions 
and treat general potentials near a minimum of the potential, see |23j . 



For t large one can use compositions and the stationary phase method. The result 
obtained is then exact since the phase in quadratic w.r.t. (x, y). 




k=l 
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4.2. Degenerate singularities. As seen in section 2 it is sufficient to study 
a micro-localized problem: 

T Z0 (E c ,h,<p) = ^Tr J e itJ t<p(t)Tp w (x,hD x )eM- l -lPh)e(Ph)dt. 

R 

Here tp G C^°(T*M. n ) is micro-locally supported near zq (cf section 2). If the 
support of tp is chosen small enough it is relatively easy to obtain a normal 
form for the phase function of the FIO approximating exp(— ^P^)S(P^). 
For the convenience of the reader we recall the contributions of equilibriums 
in the trace formula. We note S(S n_1 ) the surface of S n_1 and in the next 
two propositions it is understood that conditions to (^3) are satisfied. 

Proposition 18. If xq is a local minimum the first new wave invariant 
attached to xq are given by: 

r zo (E c ,h,ip)~h?+%- n ^ hi+^A^), 

j,ieN 2 

where the Aji are some distributions. The first new wave-invariant, attached 
to the leading coefficient is: 

frf+i:-" S( ^ J ^(vT^drj J <p(u 2 + v 2k )u n - 1 v n - 1 dudv. 

§™-i R + xR + 

Observe the simplicity of the distributional-coefficient acting on (p ob- 
tained more or less by computing a volume. Up to the spherical-mean, this 
coefficient is uniquely determined by n and k. For an unstable critical point 
the situation is more complicated and we have: 

Proposition 19. If xq is a local maximum we have : 

T Z0 (E c ,h,ip)~h% + %- n £ E hi + mog(h) m A^ m (<p). 

m=0,l j,ieN 2 

If "^fc 1 ^ ^ N, the first non-trivial new wave invariant is given by : 

J \V 2k (r))\-%dr,. 

gn-1 

At first view the result seems to be the same, but the distributions W n> k 
of Proposition [19] are respectively given by : 

(W n>k ,v) = j (C+Jtlf^- 1 + C- k \t\ n S^~\(t)dt, if n is odd, 

R 

f „i£+i_l 

(W nj k,<p) = C nk / \t\_ 2h <p(t)dt, if n is even. 

k 

The other options are given by: 
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Proposition 20. If "^fc ^ G N and n is odd then the top-order coefficients 
are: 

C^ k \og(h)h^- n ^-^ J \V 2k ( V )\-%d V J \t\^~\{t)dt. 

§™-i R 
Finally, if ^J±^ G N and n is even, and we have : 

ct k h ^- n j^ J \V2 k (vT^d v J \t\^-Mt)dt. 

§«-i R 

A careful examination of the proof shows that the last case in Proposition 
[20] is similar to the first subcase of Proposition [19] (with n odd) since C^ k = 

C~ k . But we refer to [9] for the technical details. 

Remark 21. To emphasize the consistency of these results we mention that: 

• Cn^, C^~, are non-zero universal constants depending only on n and 
k. See [9] for an analytic formulation. 

• Such terms h a and h a log(/i), a G Q never appear if E is regular. 

In this work we will mainly use the order w.r.t. h of these coefficients 
(and the constants appearing in the expansions). For a detailed proof of 
Proposition [T8l see [8] and for Propositions [191 and 1201 see [9]. The case 
k = 1, i.e., quadratic singularities, can also be retrieved from certain results 
of [1] with some support restrictions but, again, this is sufficient to attain 
our objectives. 

Asymptotic expansion at a regular energy level. With (^5) and when 
the energy E is regular, we have: 

ul-n 00 

T(E,h,<p) ~ ^LVol(S E )^(0) +J2h 1 ~ n+i c j (<p)(0) 

pg£ B i=o 

We refer to [27] for a proof. In the r.h.s. the sum concerns periodic orbits 
p of energy E and is finite since supp((^) is compact. Here S p , fi p and T p 
are resp. the action, the Maslov-index and the period of the closed orbit p 
and both Cj, D p j are differential operators of order j. If (p satisfies (^4) we 
have Cj(<^)(0) = for all j G N and for each regular value E G [Ei,E2\: 

00 

(18) T(E,h,<p) ~ ^ S ^ p, *Y, D pM(T p )V. 

peT, E j=0 

We accordingly obtain that this term is bounded and a fortiori: 

(19) T(E,h,<p)=0(l), VEe[E 1 ,E 2 ]\{E 1 c ,...,E l c }. 
This point will justify Corollary [9] 
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Next, by Lemma HH we have T p >T uniformly w.r.t. E G [E\, E?\. Hence if 
E is not critical and if in addition (H2) is satisfied the sum over the periods 
of Eq. (fT8j) is simply and in Eq. (fT9j) we obtain in fact a bound 0(h°°) (a 
fortiori (71 5) is not required in that situation). 

Asymptotics at a critical energy level. For E = E™ critical there is 
always a continuous contribution w.r.t. t in the spectral distribution showing 
up the presence of a new wave invariant. A fortiori, a choice of <p flat at the 
origin, or with a small compact support, does not erase this term. We have: 

T(E^,h, V )^J2fo(h), 
i=i 

where N m is the number of equilibrium points on S^m and each fj(h) is 
given by the leading term of Propositions 1181191 and [20l ■ 

Note that the bottom of a symmetric double well (degenerate or not) gives 
a similar answer as a single well of same nature. Hence without extra micro- 
local considerations (e.g. suitably localized eigenfunctions estimates) it is 
difficult to distinguish these 2 different settings. 

Proof of Theorem First, the micro-local Weyl-law for regular energies: 

T(E,h,tp) ~ (2 7 r/ l ) 1 - n ^(0)Lvol(S E ), 

computes the dimension n. Now assume given a critical value E c with a 
single critical point. The only choice of the spectral function tp allows to 
detect E c via the singularity f(h) of Theorem [H The knowledge of f{h) 
determines the order of the contribution. For example, if : 

f(h)~Ch a log(h), 

the critical point is a maximum and a computes the degree 2k of the singu- 
larity. With 0, the knowledge of k allows to compute the quantity : 

J \t\ nh ^- 1 ^{t)dt. 

R 

A fortiori C determines the average of |T^fc|~ on S n_1 . Without log(/i), 
the nature of the critical point can be detected by a symmetry argument 
w.r.t. tp since we a priori know n and k. In view of Propositions I18|19|20l 
we can choose tp odd, even, symmetric or non-symmetric w.r.t. the origin 
to conclude. Note that if tp is not even tp is a priori complex valued. ■ 

Remark 22. The spherical average of is a Jacobian in polar coordinates 
around xq. For example, by composition with e~'*' we obtain: 

J e -\v^)\ dx = l-T{^) f \V 2k ( V )\-^drj. 

R n S"" 1 

The same result holds by integration of /(X^fe(x)), if / G L 1 (W + ,r2k~ 1 dr). 
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Remark 23. Enlarging the list of singularities would provide a bigger "dic- 
tionary". The case of non-homogeneous singularities for V is still an open 
problem, in particular because the determination of an explicit asymptotic 
expansion w.r.t. h can be a difficult analytic problem. 

Remark 24. There are certainly many easy generalizations of Theorem [7] 
to integrable Schrodinger operators, e.g. when V(x%,X2) = V\(xi) + V^a^)- 

5. Extensions. Examples 
In this section we propose now several generalizations of the main results. 

5.1. Operators with sub-principal symbols. We will show, shortly, how 
to extend the result of Theorem [5] to the case of an /i-admissible operator. 
I mention several obvious motivations to problems involving operators with 
non-zero sub-principal symbols. 

Witten-Laplacians. B. Helffer &: J. Sjdstrand for Witten Laplacians, see 
e.g. [21j for an overview and references, have obtained recently many in- 
teresting results for these operators. For example, the Witten-Laplacian on 
zero-forms attached to the measure e~f l h is: 

A ( f °i = -h 2 A + \\Vf(x)\ 2 - ^Af(x), f e C°°(R n ), 

whose symbol p(x, £) = po(x,€) + hp\(x,£) depends on h. 

Schrodinger operators on a manifold. A Schrodinger operator at- 
tached to a Laplace-Beltrami operator on a Riemannian manifold M and 
/i-quantized by exterior multiplication: 

P h = -h 2 A M + V(x), V e C°°(M), 

generally involves a sub-principal symbol. In local coordinates with a metric 
G = G^ 1 = gij and g = det G the operator is: 

^ dxi J~q* % * dxj 

i j J J V •* J 

Hence, in the sense of the fa-calculus, we have ph := po + hpi with: 



Po(x,0 = ^2gij(x)£j£i + V{x), 




Observe that p\ = at every point where £ = 0. 
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General case. More generally, it is possible to consider /i-admissible oper- 
ators Ph whose symbols are given by asymptotic sums pt ~ ^ h 3 pj (e.g. in- 
terpreted as a Borel sum w.r.t. h) with principal symbol po(x, £) = £ 2 + V(x) 
and a subprincipal symbol p\ 7^ 0. Of course in the formula for po you can 
also replace £ 2 by the metric at x if you want to do a similar construction 
on a Riemannian-manifold (non- necessarily compact). 

Modification of the first transport equation. In the previous example 
the sub-principal symbol p\ was non-zero. This requires a light correction. 
Starting from the results of section 3 we proceed as follows. To each element 
Uh of I(K n , A) we can associate canonically a principal symbol e~h S a princ (uh), 
where 5 is a function on A such that £dx = dS on A. In fact, if Uh can lo- 
cally be represented by an oscillatory integral with amplitude a and phase 
ip, then we have S = S v = (poi^ 1 and ^^(u^) is a section of |A| 2 ®M(A), 
where: 

• M(A) is the Maslov vector-bundle of A. 

• |A| 2 is the bundle of half-densities on A. 

When p\ 7^ 0, in the global coordinates (t, y, rj) on A, the half-density of 
U h (t) is given by : 

f 

(20) v(t,y,rj) = exp(i J s (y , -ri))ds)\dtdydr}\^ . 



For this expression, related to the resolution of the first transport equation 
for the propagator, we refer to Duistermaat and Hormander [15J. Accord- 
ingly, the F.I.O. approximating the propagator has the amplitude : 

t 

d(t, z) = a(t, z) exp(i J p\(^> s (z))ds). 


Since zq is an equilibrium we have pi(<& s (zo)) = pi(zo), Vs, and : 

(21) a(t,z ) = <p(ty tpi(zo) - 

a) I£ pi(zo) = 0. This happens in many interesting situations (in particular 
for a Laplace-Beltrami operator, see above). Here the top order coefficients 
in the trace formula remains the same, also for the new wave-invariants at 
a critical energy level. 

b) If pi(zo) 7^ 0. By Fourier inversion formula we simply replace (p(t) by 
f(t + Pi(zq)) in all integral formulae of Propositions 1181191 and [201 Note 
that when using (^4), this has absolutely no effect for the mean values and 
hence on the detection of the critical energy levels. 

5.2. Eigenfunction estimates approach. We inspect now the case of an 
energy surface supporting more than one critical point. The method we use 
here is in reality much more restrictive (physically and also from the point 
of view of spectral theory) since it implicitly use eigenfunctions estimates 
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d = -Md(zi,Zj), 



via a C^°(T*M)-observable. These observable are bounded operators on 
L 2 (M) (e.g. via a Calderon-Vaillancourt estimates) and can be inserted in 
the trace, trace class operators being an ideal. 

Since everything below is local we can freely assume that M = M. n , if not we 
can use local coordinates given by the exponential map. Let K = p _1 (J) C 
T*M n be compact and: 

1 

where d is any distance on T*M n . By construction, each open ball B(z, do) C 
T*R n contains at most 1 critical point for each z G K. Clearly, we can cover 
a compact neighborhood of K by a finite number of balls B(z, do). With a 
partition of unity adapted to this covering we obtain : 

N 

Y J ^J(x,hD x ) = Id, on q?(K). 

3=1 

For each energy E 6 J, we obtain: 

Tr J 0(t)@(P h )ei t ^-^dt = J2 r & J m^^,hD x )@(P h )e^- E ^dt. 

K ■ J=1 R 

Note that the r.h.s. is studied in section 2. By the same argument as before, 
if T, s n supply) contains no critical point we obtain : 

Tr J mi>7(x,hD x )Q(P h )ei t (r h -E)dt = 0(h c °). 

And if there is exactly one critical point zq G T,e c in supp^j) we have : 
Tr f 0(t)^(x,hD x )@(P h )e^- E ^dt = ^ j (z i) )f j (h), 



and by construction no cancellation can occur since ipj(zo) > 0. 



Remark 25. In Corollary [9] we have considered (Hs) for the flow. A similar 
result holds for a chaotic dynamics and an isolated degenerate closed orbit 
can be treated as in [28J. Finally, using the results of [27] one can extend 
Corollary [9] to the case of families of periodic orbits of dimension d <n. 

It is important to notice that to put a pseudo-differential operator in the 
spectral estimates means that the inverse spectral problem is now implicitly 
expressed in terms of some L 2 -expectation: 

f j,^ j ) = (^f(x,hD x ) ( pi), 
P W \ = \ k (h)ip%, X k (h) e[E-ch,E + ch], 
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attached to eigenvectors of Ph, see, e.g., the section 'eigenvector es- 
timates' of pQ. Recall that, combining Egorov's theorem and Calderon- 
Vaillancourt estimates, in the regime h -> 0+ the measure becomes more 
and invariant under the Hamiltonian flow. It is also sometimes possible to 
obtain a measure concentrated on the critical-set (see [4]). There is no para- 
dox here since the flow is constant on the critical set of the principal symbol 
ofP h . 

For a Schrodinger operator whose potential is not a Morse function the 
inverse spectral problem seems to be ill-defined and requires eigenvectors 
estimates (which are indeed much stronger estimates than those based only 
the spectrum). Also in the previous construction an interesting problem is 
to get an a priori lower bound for the number do without doing any iteration 
on the size of supp(V^). 

5.3. Pseudo-differential operators. We can also apply the previous strat- 
egy for an /i-pseudo-differential operator with an isolated homogenous singu- 
larity as this was considered in [TO]. In fact, we will stick to the simpler case 
of a local extremum as considered in [7j. Assume that p ~ pq + hp\ + 0{h?) 
and that po has a unique critical point zq = (xo,£o) on the critical energy 
surface S^ c . Also near zq we have a conical (homogeneous) singularity: 

N 

p (z) = E c + + 0(\\(z - z )\\ N+1 ), k > 2, 

j=k 

where the functions pj are homogeneous of degree j w.r.t. z — zq. Further- 
more, assume that zq is a local extremum of pk- This implies that the first 
non-zero homogeneous component pk is even and is positive or negative def- 
inite. A fortiori zq is isolated on T,e c - An elementary example in dimension 
1 is: 

p(x,0 = ±(e + x 4 ) + R(x,0, 

where R(x,£) = 0(||(x,£)|| 5 ) and R is chosen so that p is confining and 
with tempered growth. Applying the results of [7] or |20| we can retrieve 2 
invariants of p: 

Proposition 26. The new wave invariants at E = E c determines: 
(22) A( P0 ) = ^ J \ Ph (e)\-^d6, 

§2n-l 

and the degree of homogeneity k o/p/%. 

Proof. The method of proof is the same as before, just select <p such that 
all the usual wave-invariants: 

• Energy-surface distributions, 

• Distribution supported by periodic-orbits, 
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disappear from the spectral estimates T(E,h,<p). For the new-wave invari- 
ants, following the construction of [7] or [20], it is possible to transform 
locally the phase of our FIO into —tp k - The new asymptotic problem has 
the form: 

? Z0 (E c ,h,<p) ~ j^j^ J e-^^w(t,x,Odtdxd^ w G C °°. 

RxT*R™ 

After integration w.r.t. t the becomes elliptic and a little discussion con- 
cerning the asymptotic behavior of oscillatory integrals show that: 

In 2n4-l 

T Zo (E c ,h,<p) = h~ n K 0tk {ip) + O{h—- n ). 

This result determines the even integer k. It is possible to give a full as- 
ymptotic expansion (in powers of h^) but the first new wave-invariant at 
E = E c is explicitly given by the distribution: 

1 / 2n ~ k \ 1 f _2n 

A 0ik (ip) = -(<p(t + Pl {z )),t Z() k J \Pk(0)\ ~d9, 

§2n-l 

with t zo = max(t, 0) if zq is a minimum and t zo = max(— t, 0) for a maximum. 
The lrst coefficient depends on <p, the dimension n, the degree k E 2N* 
of the singularity and the nature of the extremum. The 2nd coefficient, 
independent of tp, determines the spherical mean of p k - I 

Two different spectral-estimates, with too functions ip\ and ip 2 determine 
the ratio: 



<Pi(t + pi(zo)),t 



k 

20 



<P2(t+Pl(zo)),t 



k 

20 



so that the spherical mean on the sphere can generally be determined after 
2 spectral estimates. 

Note that many symbols would give the same value. In particular symbols 
conjugated by a rotation around the critical point cannot be distinguished 
from this new- wave invariant. 

An invariance under re-scaling. Now when zq is a minimum, we observe 
that for any j G N* we can re-scale our operator via: 

Q hJ = (p h - E c y , 

so that zq is still a minimum of the principal symbol qo = (po — E C ) J of Qh,j- 
Since: 

q0 (z) = (p k y(z) + O(\\z-z \\ k ^ 1 ), 
by a new application of Proposition 1261 at the critical value zero, we get: 



S2n-1 
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It follows that A(po) is also invariant under re-scaling (this argument does 
not apply for a local maximum) . 
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As a final remark, I mention that there is a lot of information to retrieve from 
the spectral data of certain particular Schrodinger operators on compact 
surfaces or manifolds. Certain potentials like smoothed geodesic distances 
(distance functions are singular at conjugate points) or height functions 
might describe nicely the manifold. We plan to investigate these kind of 
inverse spectral problems in a future article. 
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